Towards optimized suppression of dephasing in systems subject to pulse timing constraints 



Thomas E. Hodgson, 1 Lorenza Viola, 2 and Irene D'Amico 1 

'Department of Physics, University of York, Heslington, York, YO10 5DD, United Kingdom 
2 Department of Physics and Astronomy, 6127 Wilder Laboratory, Dartmouth College, Hanover, New Hampshire 03755, USA 

(Dated: December 10, 2009) 

We investigate the effectiveness of different dynamical decoupling protocols for storage of a single qubit in 
the presence of a purely dephasing bosonic bath, with emphasis on comparing quantum coherence preservation 
under uniform vs. non-uniform delay times between pulses. In the limit of instantaneous bit-flip pulses, this is 
accomplished by establishing a new representation of the controlled qubit evolution, where the resulting deco- 
herence behaviour is directly expressed in terms of the free evolution. Simple analytical expressions are given 
to approximate the long- and short- term coherence behaviour for both ohmic and supra-ohmic environments. 
We focus on systems with physical constraints on achievable time delays, with emphasis on pure dephasing of 
excitonic qubits in quantum dots. Our analysis shows that little advantage of high-level decoupling schemes 
based on concatenated or optimal design is to be expected if operational constraints prevent pulses to be applied 
sufficiently fast. In such constrained scenarios, we demonstrate how simple modifications of repeated periodic 
echo protocols can offer significantly improved coherence preservation in realistic parameter regimes. 
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I. INTRODUCTION 

The ability to effectively counteract decoherence processes 
in physical quantum information processing (QIP) devices is 
a fundamental prerequisite for taking advantage of the added 
power promised by quantum computation and quantum sim- 
ulation as compared to purely classical methods. Dynami- 
cal decoupling (DD) techniques for open quantum systems 12 
have been shown to be able to significantly suppress non- 
Markovian decoherence for storage times that can be very 
long relative to the typical time scales associated with the 
decoherence process itself. Over the last decade, the de- 
sign and characterization of viable DD schemes for realis- 
tic qubit devices has spurred an intense theoretical and ex- 
perimental effort, taking DD well beyond the original nu- 
clear magnetic resonance (NMR) setting 1 . While earlier DD 
schemes relied on the simple periodic repetition of instan- 
taneous pulses (so-called 'bang-bang' periodic DD, PDD 2 , 
and its closely-related time-symmetrized version, so-called 
Carr-Purcell DD, CPDD 3 ^ 4 .), recent theoretical investigations 
have explored the benefits of more sophisticated control de- 
sign in a number of ways. In particular, this has led to de- 
vising recursive and randomized pulse sequences for generic 
decoherence models on finite-dimensional quantum systems 
- so-called 'concatenated' DD (CDD 5 ^ 6 -) and 'randomized' 
DDiS; to identifying 'optimal' protocols for a single qubit 
undergoing pure dephasing - most notably, the so-called 
Uhrig DD (UDD) 9 - 10 ' 11 ' 12 ! 13 ' 14 , and its extension to 'locally 
optimized'^i£ DD sequences tailored to specific noise envi- 
ronments; and, most recently, to combining the advantages 
of concatenation and optimization for a single qubit exposed 
to arbitrary decoherenc o 17 ' 18, 19 . As a key common feature, 
these investigations highlight the sensitivity of DD perfor- 
mance to the details of the applied control path, and point to 
the importance of carefully tuning the relative pulse delays 
in order to boost the efficiency of the achievable decoherence 
suppression 2 ^. 

In view of the above rich scenario, assessing the perfor- 



mance of different DD protocols in specific qubit devices 
and/or in the presence of specific control constraints becomes 
especially important. Recently, the effectiveness of traditional 
multi-pulse spin-echo sequences based on PDD and CPDD, as 
compared to 'high-level' protocols based on CDD and UDD, 
has been scrutinized in several control settings. In particu- 
lar, a number of theoretical studies have addressed suppres- 
sion of pure dephasing associated to spectral diffusion^ 3 - and 
hyperfine-induced decoherence^ 4 . from a quantum spin-bath 
for an electron-spin qubit, as well as suppression of classical 
1/ / phase noise in a superconducting qubi t 25,26 . Experimen- 
tally, the performance of CDD protocols has been character- 
ized for an NMR spin qubit 27 , while optimal UDD implemen- 
tations have been reported for both a trapped ion qubit ex- 
posed to engineered classical phase nois o 15,16,28 and, most re- 
cently, for electron spin qubits undergoing spin-bath decoher- 
ence in a malonic acid crystal 2 ^. These studies have demon- 
strated, in particular, how UDD can significantly outperform 
low-level DD schemes provided that the noise spectrum has a 
sharp high-frequency cutoff and sufficiently high pulse repe- 
tition rates may be afforded. 

Amongst prospective solid-state QIP platforms, exciton 
qubits in self-assembled quantum dots (QDs) have likewise 
received vast attention in recent year a 30 ' 31 : due to the cou- 
pling to photons, excitons can be driven ail-optically on sub- 
picosecond time scales 30 . Excitonic implementations also al- 
low the flexibility of designing hybrid solid state-flying qubit 
scheme a 32 ' 33 . Pure dephasing turns out to be the dominant 
factor limiting the coherence lifetime in such qubit devices, 
where strong coupling with phonon modes of the host crystal 
result in typical decoherence (T2) time scales of a few pico- 
seconds^ 4 .. We have previously shown in Ref. [HI that, remark- 
ably, PDD allows for substantial exciton coherence recovery 
in experimentally relevant parameter regimes (up to 90% re- 
covery over ~ 10 ps at room temperature), the control per- 
formance being especially enhanced for QD shapes and bias 
fields optimized for quantum computing architectures. Our 
goal in this paper is to quantitatively assess to what extent 
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more elaborated DD schemes - in particular, sequences em- 
ploying non-uniform pulse timings - can improve beyond the 
simplest PDD setting when a lower bound on the achievable 
control time scale (minimum pulse separation) is present. 

We find that in the presence of such a timing limitation, 
simple protocols such as PDD or CPDD may outperform high- 
level sequences based on CDD/UDD. Interestingly, on the one 
hand this reinforces similar conclusions drawn in Ref. |26| for 
classical dephasing in superconducting qubits. On the other 
hand, we additionally show how it is possible to engineer 
a suitable 'preparatory' sequence that enhances the perfor- 
mance of a subsequent PDD pulse train. In the process, we 
take advantage of the exact solvability of a purely dephas- 
ing model in the presence of instantaneous pulses to obtain 
an exact representation of the controlled dynamics in terms of 
the free evolution. This allows rigorous results on the long- 
time asymptotic decoherence behaviour to be established for 
generic noise spectral densities, by allowing in particular a 
comparison between ohmic and supra-ohmic environments. 
Furthermore, our work provides a first explicit analysis of 
CDD performance in the presence of a quantum bosonic bath. 
From a practical standpoint, our results suggest that simple 
DD protocols may remain a method of choice if significant 
timing constraints are in place, and that incorporating such 
constraints from the outset is necessary before further opti- 
mization can show its benefits. While our numerical results 
are tailored to excitons in QDs, we expect the above conclu- 
sions to be relevant for other constrained qubit devices. 



II. SINGLE-QUBIT DEPHASING DYNAMICS 

We consider the pure dephasing dynamics of a single qubit 
coupled to a non-interacting bath of harmonic oscillators. The 
Hamiltonian of such a system may be written in the form 

H = —a z + h^Wjtfbj 

3 

+ hY,(9p] + 9jbj)[(l - a)a Q + aa z ] (1) 

3 

= H + hJ2(9j b ]+9jbj)[(l-a)<ro+aa z ], (2) 

3 

where E gives the energy difference between the qubit's lev- 
els, bj and bj are canonical creation and annihilation opera- 
tors of the oscillator mode j, and gj describes the coupling 
between the qubit and the j-th bath mode. In the above ex- 
pression for H, the parameter a accounts for the possibility 
that either both or only one of the spin (or pseudo-spin) qubit 
computational levels effectively couple to the bath: a = 1 cor- 
responds to the standard purely-dephasing spin-boson model, 
whereas if a = 1/2, only the a z = +1 eigenstate couples to 
the bath. Specifically, for an excitonic qubit, the logical states 
are represented by the presence or absence a single (ground- 
state) exciton in the QD^2, and E is the energy relative to the 
crystal ground state. 

As time evolves, the qubit becomes entangled with the en- 



vironment and the off-diagonal elements of the qubit density 
matrix evaluated at time t in the interaction picture with re- 
spect to Hq readied 

poi(i) =Pw(t) =poi(* = 0)e- r W (3) 

r(t) = r (t) (4) 
, N9 f°° , i(u) , / hv \ r , N1 

= (2a) 2 ^ d^coth(— )[l-cosM)], 

where T is the temperature, ks the Boltzmann's constant, and 

I(u,) = Y,S(u;-u )j )\g j \ 2 (5) 

3 

is the spectral density function characterizing the interaction 
of the qubit with the oscillator bath. For a supra-ohmic envi- 
ronment, I(uj) a; 3 , as opposed, for instance, to an ohmic 

reservoir where I(uj) w ~ u>. Likewise, the high-frequency 
behaviour is characterized by a frequency cut-off oj c , for in- 
stance, for excitons one can assume that I(uj) e~ w . 

As it turns out, the decoherence of the qubit in the pres- 
ence of an arbitrary sequence of bang-bang pulses, each ef- 
fecting an instantaneous 7r rotation, can still be exactly de- 
scribed by Eq. (f3]), provided a modified decoherence func- 
tion is use d 2, 10,36 . Consider an arbitrary storage time t, dur- 
ing which a total number s of pulses is applied, at instants 
{ti, ...,£„,..., t s }, with < ti < t 2 < ■ ■ ■ t s < t. By using 
the theory developed by Uhrig in Refs. [TolfTll we can define a 
controlled coherence function T(t) in the following way: 

f r (t) t < h, 
r(t) = <^ r n (t) t n < t < t n+1 , o < n < s , (6) 
[ r,(t) t s < t. 

Here, To(t) is given in Eq. (0]i whereas for 1 < n < s we 
le02 

r„(t) = (2«) 2 ^ °^coth(^)| y „M)| 2 ^, 

n 

y n (z) = 1 + (-l)" +1 e" + 2 ]T (-l) m e" s "% z > 0, 

m— 1 

with the n-th pulse being understood to occur at time t n = 
8 n t, and < Si < . . . S n < . . . S s < 1. While the instanta- 
neous pulse assumption must be handled with care in general, 
we have discussed in Ref. [35] how it translates into reason- 
able physical constraints for an excitonic qubit coupled to a 
phononic bath. 

We now proceed to directly relate T n (t) to To(t) for ar- 
bitrary n. Let us first rewrite the above coherence function 
T n (t) in a compact way as 

/>oo 

r„(t)= / 7 ? (^)|?;„M)| 2 dw, n>0, (7) 
Jo 
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where we have defined 

|y M)| 2 ^|l-e^| 2 , 

and 

By relating \yi(cut)\ 2 to \yo(u>t)\ 2 we can write 

rx(t) = -r (t) + 2r (ti) + 2r (i - t x ). 

Upon continuing this iteration we find 

r 2 (t) = -ri(t) + 2ri(t 2 ) + 2r (t-ti) I 
r»(t) = -r B _i(t) + 2r n _i(t„) + 2r (t-t„). (ii) 

Furthermore, by expressing |y„(u;i)| 2 as a function of 
|yo(wt) | 2 , we are able to write the entire evolution in the pres- 
ence of an arbitrary pulse sequence only in terms of the un- 
controlled evolution. Explicitly, we find: 

n 

r„(t) - 2^(-i) m+1 r (t m ) 

m—1 

n 

+ 4 E E r o(*»-*i)(- 1 ) m " 1+J 

m— 2 j<m 

+ 2 ^ (-l) m+n r (< - t m ) + (-l)'To(i). (12) 

m— 1 

The above equation is one of the main results of this paper. By 
using Eq. (1I2L it is, in particular, straightforward to see that 

Tn-^tn) = Hm r n (t) = r„(t„). (13) 

This confirms that the function T(t) as defined in Eqs. (O is 
continuos at the (instantaneous) pulse timings, as expected on 
physical grounds. 

As a first example of the usefulness of this representation, 
we consider how two pulses may be used to increase the 
asymptotic coherence of a supra-ohmic system, in which the 
free dephasing dynamics saturates in the long-time limit to 
a finite valu o 34 ' 37 ro(oo) > 0. Taking the t — > oo limit in 
Eq. ( [Tot or, equivalently, letting n — 1 in Eq. (fT2l . yields 
ri(oo) = 2r (<i) + r (oo). Since T (t) > for all t, this 
shows how a single pulse cannot decrease the asymptotic de- 
coherence level. However, after two pulses we have 

r 3 (t) = ro(t)-2r (t-ti)-2ro(t 2 ) 

+ 2r (t 1 ) + 4r (t 2 -ti) + 2ro(t-t 2 ). (14) 

Therefore, 

r 2 (oo) = r (oo) - 2r (t 2 ) + 2r (ti) + 4r (i 2 - t x ), (15) 

and ti and i 2 can be chosen to decrease the asymptotic deco- 



FIG. 1: Comparison between Fo(t), ^i(t), and Ti(t) for an exciton 
qubit at T = 77 K, as computed from Eq. ((6). Pulse times are 
t\ = 0.2 ps and = 0.31 ps. 



herence provided that 

2r (t 2 )-2ro(ti)>4ro(*2-*i). (16) 

In Fig. Q] we plot T(t) for an exciton qubit coupled to phonon 
modes and subject to two control pulses at t\ = 0.2 ps and 
<2 = 0.31 ps. For comparison, we also plot the evolution 
under a single control pulse at t\ = 0.2 ps and the free evo- 
lution r (t). As one can see, Eq. (fToT ) can indeed be satis- 
fied. Numerical results showing how a few pulses can increase 
the asymptotic coherence have been reported for excitonic de- 
phasing in Ref. [3H 

For the case of Fig. Q] as well as for all the numerical ex- 
amples in this paper, we consider (unless otherwise stated) 
an exciton qubit tightly confined within a GaAs QD at 77 
K. The QD potentials are modeled as parabolic in all three 
dimensions, with confinement energies in the z-direction of 
hid e = 505 meV and hwh = 100 meV, while hio e = 30 meV 
and hwh = 24 meV in the in-plane direction a 30 ' 39 . The sub- 
script e/h indicates electron/hole, respectively. For this exci- 
ton, in the absence of control most of the coherence is lost af- 
ter a few picoseconds 39 . Having this specific system in mind, 
we shall set a = 1/2 henceforth in our numerical calcula- 
tions, and plot the quantity | exp(— T n (t)) | 2 , which is directly 
proportional to the square modulus of the measured optical 
polarization P(i). 

As discussed in detail in Ref. [35L the spectral density of this 
system is given by 

Hu)=I e (u) + I h (u) + I eh (u), (17) 

where the indices e/h/ 'eh correspond to single particle spec- 
tral densities of the electron and the hole, and to the electron- 
hole inteference term respectively, and 

h/h/M = E n' hleh ^) exp ( - a"' ) ■ (18) 

j V Ci.e/h/eh ' 

Here, i labels different phonon modes, whereas F^' h ' eh (uj) 
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is a mode-dependent function for which p^l h l eh ( w ) < 
The spectral density may be further approximated as 

2 

I(lu) «F W 3 exp(-^), (19) 

where the parameters F and uj c are determined by fitting a 
curve of the form Eq. ( fT9l ) to the actual exciton spectral den- 
sity. For the particular exciton parameters listed above, this 
yields F = 1.14 x 10 _26 s and hu c = 2meV. 



III. PERIODIC DD: PERFORMANCE AND EXACT 
ASYMPTOTIC PROPERTIES 

For a Hamiltonian as in Eq. (Q]), a DD cycle consisting of 
two uniformly spaced rotations by tt about the x axis, 

XAtXAt, (20) 

where time ordering is understood from right to left, removes 
the interaction between the qubit and the boson bathed to the 
lowest (perturbative) order in oj c T c , with T c = 2At. The sim- 
plest DD protocol, PDD, is obtained by iterating the above 
control cycle in time. 

Fig.|2]compares the free evolution with the PDD-controlled 
dephasing for the exciton qubit under examination, computed 
from the exact expressions given in Sec [TT] Sequences with 
three different pulse delays are shown, At = 0.1 ps, At = 
0.2 ps, and At = 0.3 ps, respectively. For the exciton qubit, 
two conditions determine a suitable range of At for effective 
PDD: i) On the one hand, it is necessary that the control time 
scale T c be sufficiently short with respect to the (shortest) cor- 
relation time of the decoherence dynamics, which means in 
this case 2 At < r c = 2ir/u c . Physically, this can also be 
interpreted by requiring that the characteristic frequency in- 
troduced by the periodic control, 

Wres = A? 

be significantly higher than the spectral cut-off fequency it- 
self, w res > u>o m sucn a way that the DD-renormalized 
spectral density function, I(w) tan 2 (to At/2), is effectively 
'up-shifted' beyond the bath cutof f 12 i 19 i 35 i 41 . ii) On the other 
hand, the existence of a lower bound on the pulse duration 
implies a lower bound on the separation At in order for the 
instantaneous-pulse description to be accurate. As discussed 
in Ref. [35[ this means At > 0.1 ps for semiconductor self- 
assembled QDs of interest for QIP 

The values of At used in Fig. [21 are consistent with both 
these conditions. It can be seen that coherence decays until the 
first bit-flip occurs, after which it rises, reaches a local maxi- 
mum before decohering once again - with this pattern repeat- 
ing between every two bit-flips. It can also be seen that DD 
recovers most of the dephasing, that is, exp(— T(t)) is much 
closer to unity than in the uncontrolled evolution, which falls 
rapidly before saturating to cxp(— r (oo)). After the first few 
initial pulses, the dephasing enters a phase in which the be- 
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FIG. 2: | exp(— F(t)) | 2 for the exciton qubit in the presence of PDD 
with At = 0.1 ps, At = 0.2 ps, At = 0.3 ps, compared with the 
free evolution determined by To(t). 



haviour after the (n + l)-th pulse is approximately the same as 
the one after the n-th pulse. For At = 0.1 ps and At = 0.2 ps, 
the average dephasing over each cycle in this 'steady-state' 
phase is very small, leading to a practical 'freezing' of the 
average decoherence over a period much longer than the es- 
timated (sub-picosecond) gating times 3 ^. For At = 0.3 ps, 
however, the increase of decoherence due to this average de- 
phasing with time is more noticeable, leading to worse DD 
performance overall. It can also be seen that, to minimize the 
effects of dephasing, any readout on the qubit should be made 
half-way between two control pulses. As it is well known in 
NMR, this motivates a proper choice of the observation win- 
dow, which underlies the Carr-Purcell (CP) sequenced and is 
also discussed in Ref. |42j in the spin-boson context. 



A. Long-time dynamics: Ohmic versus supraohmic behaviour 

A main advantage of the exact representation established 
in Eq. ( flZb is that it allows detailed quantitative insight on 
the controlled dephasing behaviour to be gained. In particu- 
lar, we focus on long-time coherence properties, which have 
also received recent attention in view of control-dependent 
'saturation' effects observed in the context of spin-bath 
decoherence^ 3 - (see also Ref. l44b . We start by quantifying 
how the decoherence function in the presence of n pulses dif- 
fers between two consecutive control times. Let 

Ar„ = r„(t„ +1 ) -r n _i(t n ). (21) 

By using Eq. ( fT2l we obtain: 

Ar„ = (-i)"[r (t„ +1 )-r (t„)] + 

n 

+ 2^r (t„ +1 -t i )(-i)" +j 

n-1 

- 2^r (t n -t i )(-i) J+n - (22) 
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FIG. 3: Differential dephasing function, Ar^ DD , for the exciton qubit 
under examination in the presence of PDD with At = 0.3 ps (top) 
and At — 0.25 ps (bottom), calculated from Eq. d22t . The dot- 
ted lines show, in each case, the limiting value AT^ given by Eq. 
l |23| >. Notice that for n < n s a, where n sat ~ 15, the sign of Arf t DD 
oscillates, in agreement with Eqs. iB5\ and dB7t . 



Let now Arf t DD denote the above 'differential dephasing 
function,' Eq. (12 It . specialized to a PDD protocol. Then, as 
showed in AppendixlAl the following asymptotic result holds 

for an arbitrary dephasing environment: 



AT 



lim AT; 

n— ►oo 



8w res r? (w res ) . 



(23) 



Interestingly, Eq. (|23l can be used to describe how the de- 
phasing function changes between any two instants separated 
by At, for large enough t. That is, consider 

Ait DD (f) = r n+1 (t + t Il+1 ) - r„(t + t n ), (24) 

where < i < At, t n = nAt. By using Eq. (fT3l we can 
verify that ArP DD (0) = Ar™ D . Then one may also prove 
(see Appendix [B] for detail) that 



Arr (f) n ^> 



AT r 



(25) 



where n sat = t sat /At is a sufficiently large integer defined in 
the same Appendix. Eq. (l25T l shows that the dephasing in- 
crement becomes independent of n and t for t > t sat , that 
is, dephasing asymptotically enters a periodic oscillation 'in 
phase' with the PDD sequence. Thus, Ar^ in Eq. (|23T > may 
be used to describe the difference in dephasing between any 
two times separated by At - in particular, between consecu- 
tive coherence maxima which for t > t sat occur at t « At/2. 
For a supra-ohmic environment as in the exciton qubit, the 
convergence of Ar™ D to Ar^, Eq. (|23l , is very fast. This is 
illustrated in Fig. [3]fbr two representative values of At. 

Because Ar,^ in Eq. d23l is non-zero as long as At is fi- 
nite, we can infer that r n diverges for fixed At as n — > oo. 



While this in principle implies a decay of exp(— T(i)) to zero 
under the PDD, details of the spectral density function (in- 
cluding the nature of the coupling spectrum and the form of 
spectral cutoff) become essential to characterize different dy- 
namical regimes of interest. In what follows, we illustrate 
these features by contrasting ohmic and supraohmic dephas- 
ing environments, and by considering stroboscopic sampling, 
t n = 2nAt, in which case explicit analytic expressions for the 
PDD 'filter function' \y2 n (2nujAt)\ 2 are available. Specifi- 
cally, upon combining Eq. (1 lb) of Ref. 07J with Eq. ( [TBI re- 
covers the well-known resul t 1,1 1,42 : 



T 2n (2nAt) 



9 fUjAt\ 

477(0;) sin (wnAt) tan ' 



2 r- 



(26) 



In general, we expect two dominant contributions to the 
above integral: the one from small values of uj, where rj(u>) 
is not small, and the one from the region of the resonance, 
uj ~ u; re s, where |?/2n(wt)| 2 may be large. First, note that for 
both a ohmic and supra-ohmic spectral density, the contribu- 
tions from the small-w region saturate to a. finite value with 
time. For the ohmic case, this is true irrespective of the fact 
that the free dephasing dynamics does not exhibit a similar 
long-time saturation. This behavior is due to the control term 
tan 2 (u At/2), which increases the rate at which the integrand 
goes to zero as uj — > 0. Second, the contribution from the 
<jj w oj ks region is more or less relevant depending on the form 
of the spectral cutoff. Clearly, such 'resonating' contributions 
do not pose a problem in the limiting situation of an arbitrarily 
'hard' spectral cutoff of the form O(o; — ui c ) (9( ) denoting 
the step function), since, as remarked earlier, o; res > lo c in a 
good DD limit. For a smooth ('soft') spectral cutoff, the res- 
onating contribution increases with time and will ultimately 
be responsible for the divergence of r2„(2nAi) as n — > 00. 
In fact, AFqo corresponds precisely to such a frequency range. 
As shown by Eq. d25l l, we can approximate Ar„ w AT^ for 
t > t sa t: since at such long times, the contributions to Eq. (|26| | 
from small o> have saturated, dephasing is indeed dominated 
from the region around o> res . Thus, for both ohmic and supra- 
ohmic systems under PDD, the coherence will eventually de- 
cay to zero for large enough times and soft cutoffs. 



The above considerations are illustrated in Fig.|4] where we 
plot exact results calculated from Eq. d26l i for a representative 
ohmic spectral density with an exponential cutoff 1 : 



I a (uj) = Fuj exp 



UJ 

UJr 



(27) 



In order to highlight the different contributions to the overall 
dephasing function, we also explicitly compute and plot the 
following quantities: (i) (dotted line) 

r smtJ (2nAt) = / V (uj)\y 2n (Lo2nAt)\ 2 duj, (28) 
Jo 
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co c t 

FIG. 4: Dephasing behavior for an ohmic spectral density with expo- 
nential cutoff as in Eq. (27), with F = 0.5, a = 1/2, At = 0.0015, 
uj c = 100 and T = 100o; c , in units where H = ks = 1. While stro- 
boscopic sampling is implied, continuous interpolating lines are used 
for clarity, a) Full decoherence function, exp(— T(t„)), Eq. \26\ 
(solid line); low-frequency contribution, cxp(— T smu} (t)), Eq. {2Sj 
(dotted line); resonanting contribution, exp(— r res (t)), Eq. \29\ 
(dashed line) versus rescaled time u c t. b) Comparison between 
exp(-r re8 (2nAt)) (Eq. <(29j) (points) and expi-AF^t/ At) (solid 
line). 



which isolates the small-w contributions, and (ii) (dashed line) 



T res (2nAt) 



3oj„ s /2 



<W2 



Ti(w)\y 2n (w2nAt)\ 2 duj, (29) 



which isolates the contributions from the ui ~ w res region. 
Three distinct regions may be identified: an initial drop in co- 
herence due to the low-frequency modes, until saturation of 
Eq. d28l ) occurs at about t = r c ; a plateaux region where the 
contributions from Eq. d29l ) are not important enough to cause 
further decoherence; and a final decay of coherence to zero 
caused by increasing contributions from the u> w uj ks region. 
Fig- S] also compares (bottom panel) the resonating contribu- 
tions calculated from Eq. ([29b with the asymptotic prediction 
cxp(-Ar oo i/At) (solid line), with AI^ = 4.507 x 10~ 7 . 
The data confirm that Ar^ does indeed arise from the res- 
onating contributions as expected, and that as long as the low- 
frequency contributions have saturated, Ar^ may be used 
to accurately describe dephasing under PDD in the long time 
limit, that is, Ar„ w AT^, for t > i sat . 

Additional insight may be gained by examining how the 
above different regimes (initial decay, plateaux, final coher- 
ence decay) are affected by the harder or softer spectral cutoff 
function. Beside the ohmic spectral density of Eq. d27l i. con- 
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FIG. 5: Dephasing behavior for supraohmic spectral densities with 
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sider the following supra-ohmic spectral densities: 

I b (uj) = Fw 3 exp ( - — ), 

2 

I c {lo) = Fw 3 exp (- ^a), 



(30) 
(31) 



where, in particular, I c (oj) has a Gaussian tail, similar to the 
excitonic qubit case. When comparing /&(cj) and I c (ui) (see 
Fig. 0, the harder cut-off due to the Gaussian tail strongly 
reduces the value of r?(w res ), and hence greatly increases the 
duration of the plateaux regime. In fact, for the set of param- 
eter chosen, our numerics loose the necessary precision well 
before the third regime sets on for I c {w). The harder cut-off of 
the Gaussian case also decreases r smw and, in turn, decreases 
the decoherence that occurs before the plateux. 



B. Short time dynamics 

In the previous section we analyzed the dephasing dynam- 
ics in the presence of PDD for t > t sat . Here, we focus on 
t < t sat . The long time regime is entered when V n +i = 
r„ + AToo, and for this to occur the coherence must oscillate 
in phase with the DD pulses. However, the natural response of 
the coherence after the first PDD pulse is instead to oscillate 
with a period of 2 At (twice that of PDD pulses, recall Fig. [6]). 
This follows from the fact that the first bit-flip occurs an in- 
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terval At after a maximum, ro(0), and for sufficiently small 
At, the dephasing function is roughly symmetrical about the 
control pulse, so the coherence maximum following the first 
pulse occurs at t w 2At. The PDD sequence quickly drives 
the coherence into phase with it, (see Fig. [6), but the first few 
even bit-flips in PDD occur near the coherence maxima, and 
this worsens the performance of the control sequence. This 
may be seen by considering Eq. ( fTTT i at time t = t n + t, with 
< t < At. By expanding the first and last terms to first 
order in t, and considering that r (0) is a maximum, Eq. (fTTT i 
rewrites as 

at 

The second term in the above equation is a constant, hence 
there can be a coherence peak after the n-th control pulses 
only if the derivative dT n -i{t n ) / dt > 0, as also pointed out 
in Ref.HH In particular, again using Eq. (fTTT i. we can calculate 

dT n (tn) _ r ra ft„ + t) - T n (t„) _ dT n -l(t n ) 

dt ~ t dt 

which shows that the larger the gradient of r„_i(t„), the 
faster the coherence is retrieved immediately following the 
?ith pulse. In particular, if r n _i(i) is locally flat at the time of 
the n-th pulse, no coherence gain can occur after that pulse. 

We can see from Fig. [6] that as PDD drives the coherence 
oscillations into phase with it, AT n has alternating sign for 
odd and even n (cf. Eqs. ( IB5b and (IB7t ). Ar n is initially 
negative for odd n and positive for even n, while its magnitude 
decreases until a time t av after which Ar„ becomes positive 
for odd n and negative for even n, before saturating to AT ra = 
Ar^. We see numerically that t av is independent of At, with 
i av ~ 0.5 ps in our case. Furthermore, we can show from 
Eq. ( I A It that if we consider the times at which the control 
pulses occur (t = 0), then 

Ar™ D (0) = Ar^O) + (-l)"Ai 2 r£(nAt), (33) 

where 

d 2 T (nAt) _ r ((n - l)Ai) - 2r (nAt) + T ((n + l)At) 
dt 2 " At 2 

From this expression we can understand the behaviour of the 
dephasing for PDD as the coherence oscillations are driven 
into phase with the PDD pulses. As n increases, the sign 
of the last term in Eq. ( f33l ) alternates, and its magnitude 
decreases as dTo(nAt)/dt reaches a maximum before de- 
creasing and tending to zero (recall the behavior of To(t) in 
Fig-ID- Thus, we can now rigorously define t av by the condi- 
tion d 2 To(t av )/dt 2 = 0, that is, when the gradient of To(t) is 
maximum. 



C. Practical considerations 

Even if the qubit coherence eventually decays to zero un- 
der PDD in our excitonic system, for practical purposes we 
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FIG. 6: Short term dephasing of the exciton qubit under PDD, with 
At = 0.1 ps. The diamonds indicate the timing of the PDD pulses. 



only need to suppress the dephasing for the qubit lifetime, 
T\. From the above discussion, we can estimate more pre- 
cisely how short At must be, in order for this to happen. For 
t = nAt + At/2 > t sat , we can approximate the off-diagonal 
density matrix element at the maxima of coherence (where a 
measurement would be made) as 

p i[(n + l/2)Ai] « p 01 (o) e - r -.[( n -+ 1 /2)At]-(™-n sill )Ar M ^ 

(34) 

Considering the long-time limit, if At is sufficiently small 
and n ^> n sa t, we may neglect the coherence that is lost whilst 
t < t sat , and further approximate the dephasing as 

Poi[(n+l/2)At] «poi(0)e- nAr ~ -poiWe-^'. (35) 

Thus, in the long-time limit, we effectively have 1/T| ff = 
ATao/At. A sufficient condition for the dephasing to be sup- 
pressed for the entire qubit lifetime is then 

Tf=-^>n. (36) 

Fig. [7] shows AToo/At as a function of At for the exciton 
qubit under consideration, for which 1/Xi = 1 ns _1 . It can 
be seen that for At < 0.2 ps, PDD effectively suppresses 
dephasing for the entire lifetime. This is in excellent agree- 
ment with our previous results in Ref. HH where we found 
numerically that At = 0.2 ps leads to efficient PDD, but, in 
comparison, At = 0.3 ps could only suppress the dephasing 
for relatively short times. 



IV. COMPARISON OF PDD WITH NON-UNIFORM DD 
SCHEMES 

Having characterized the performance of the simplest DD 
scheme, where the control involves a single time scale At, 
we proceed to examine some of the high-level protocols men- 
tioned in the introduction, which involve non-uniform pulse 
delays to a lesser or greater extent. While CPDD is both, his- 
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FIG. 7: Effective long-time coherence decay rate, 1/Tf ff = 
AToo/ At, Eq. ( 1361 ), as a function of At. The dashed line is the 
qubit inverse lifetime, 1/Ti = 1 ns _1 . 



Sequence 


Pulse Timing 


So 


Free(At) 


Si 


XAtXAt 


S-2 


X[XAtXAt]X[XAtXAt] = AtXAtAtXAt 


S 3 


X[AtXAtAtXAt]X[AtXAtAtXAt] 


Se 


XSi-jXSe-i 



TABLE I: Concatenated pulses sequences for a purely dephasing 
single-qubit interaction. Time ordering is from right to left. 



torically, the most established approach and, ultimately, one of 
the most effective, we defer its discussion until after the anal- 
ysis of CDD and UDD, since it turns out that for the supra- 
ohmic system at hand CPDD naturally suggests the optimiza- 
tion strategy that will be introduced in Sec. [V] 



A. Concatenated decoupling 

Instead of repeating the basic control cycle given in 
Eq. ( |20t , CDD recursively concatenates it within itself. Let 
Se denote the sequence corresponding to the £-th level of con- 
catenation, as given in TableQ] 

For a qubit undergoing arbitrary decoherence, CDD with 
a 'universal decoupling' cycle given, for instance, by 
AtXAtZAtXAtZ, has been shown 5 to significantly outper- 
form PDD in the limit At — > 0. However, for purely de- 
phasing systems for which At has a finite lower limit, and 
for single-axis protocols constructed out of the basic cycle in 
Eq. d20l) . the advantages of CDD are largely lost, and PDD 
may be more efficient 4 '-. While different ways for comparing 
different DD protocols can be considere d 5 ' 24 ' 26 , we shall fo- 
cus here on comparing the efficiency of PDD and CDD at en- 
suring dephasing-protected storage of the exciton qubit for a 
fixed time T stora ge- In particular, for our calculations we choose 
^storage = 10 ps. This time is appropriate given the typical gat- 
ing time for exciton-based QIP, which is of the order of 1 ps^2. 



1. Single CDD cycle 

Given T storage and the presence of a physical constraint on 
At, a first way to exploit CDD is to identify a minimum con- 
catenation level, £*, for which the length of the corresponding 
sequence, Tg* = 2 e At, exceeds r storage . For a given At, 
increasing I beyond this point would not modify the results 
because the pulse timings over T stora g e would be unchanged, 
(see Table Figure [8] compares CDD and PDD for storage 
of an exciton qubit for different At. As expected from the 
general analysis of Ref. HI the efficiency of CDD increases 
with decreasing At. However, in the range of values under 
exploration, and with readout effected at the maxima of the 
coherence curve, CDD is found to be more efficient than PDD 
only if At < 0.036 ps. The latter time scale is substantially 
smaller than physically allowed in our system. 

We can understand the possible advantage of CDD by com- 
paring it with the long- and short-time behaviour of PDD 
(Sees. MI Al and I11IBI respectively). Eq. d34l ) shows that 
the long-time performance of the protocol depends on Ar^, 
and r„ sal [(n sat + 1/2) At]. For very small At (hence small 
AToo), PDD is not the most efficient scheme because it leads 
to a value of r„ sM [(rt sat + 1/2) At] which may be greater 
than for other pulse sequences, due to the initially 'out of 
phase' pulses. In the regime where CDD outperforms PDD 
(very small At), the contributions to dephasing from around 
oj = uj ks (see Sec. MI Al ) are negligible for both sequences 
over T storage , since for t > t sat both sequences preserve the 
maxima of coherence very close to the value cxp(— T(t™ x )) 
corresponding to the time i™ x of the first maximum that fol- 
lows t sat . The advantage of CDD (if any) comes from the 
different behaviour of the dephasing over the first few control 
pulses, that is, up to t = t sat . The timing of the pulses in the 
CDD sequence are similar to those of PDD, but with fewer 
pulses at the instants where the even pulses occur in PDD. 
These 'missing' pulses are those which would occur near the 
maxima of coherence in the initial stages of the sequence (see 
insets of Fig. [8]), that is, the ones responsible for decreasing 
the coherence maxima while t < t sat in PDD (Sec. IIII BD . 
These 'missing' pulses allow the dephasing to maintain its 
natural response frequency after the first bit-flip, and no loss 
of dephasing is needed to change the rate of the oscillations 
of coherence. Therefore, r CDD (t™ x ) > r PDD (t™ x ), and for 
t < T storage , T(t) w r(t™ x ) for both PDD and CDD in the 
limit of sufficiently small At. 

While the above explains why CDD may outperform PDD, 
as soon as At is long enough such that AT^ is significant 
over T stora ge, PDD becomes the most efficient sequence. The 
period of the coherence oscillations for CDD is twice the one 
for the PDD sequence corresponding to the same At (see in- 
sets in Fig. [8]), resulting in faster dephasing at long times t for 
CDD. 



2. Periodic repetition of CDD cycles 

A different use of CDD consists in truncating concatenation 
at a fixed level and periodically repeating the resulting 'super- 
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FIG. 8: (Color online) | exp(-r(t))| 2 for PDD (black line) com- 
pared with CDD (light line) for At = 0.016 ps (t = 10, top), 
At = 0.036 ps (t = 9, middle), and At = 0.055 ps (t = 8, 
bottom). Insets: Close-ups of the same evolutions at short times; 
the pulse timings are indicated as well, with crosses (CDD) and dia- 
monds (PDD). 
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FIG. 9: (Color online) Comparison of PDD (dark line) and PCDD 
(light line) protocols with At = 0.1 ps. Top: Second-level con- 
catenated cycle, PCDD2. Bottom: Third-level concatenated cycle, 
PCDD3. Inset: Zoom over the initial part of the time window with 
timings of the pulse sequences explicitly indicated (diamonds for 
PDD and crosses for PCDD 3 ). 



PDD for 1 = 2, but worse for I = 3. The difference between 
PCDD2 and PCDD3 may be understood as a consequence of 
the fact that in terms of a Magnus expansion 2 ', concatenated 
cycles with even I are time-symmetric, thus cancel the inter- 
action with the phonon bath up to (at least) the second order. 
Over the time period shown, PCDD 2 also outperforms stan- 
dard PDD (see Fig. |9j upper panel). However, the coherence 
oscillations for PCDD2 occur over a period of 2Ai since, af- 
ter the initial pulse, the sequence is equivalent to PDD with a 
base time interval of 2At. Therefore, we expect PDD to be 
more efficient for long storage times, as PCDD2 will yield a 
larger Ar^ than a PDD sequence characterized by At, hence 
worse asymptotic performance. 



cycle', constructed from Table. Q] For instance, truncation at 
1 = 2 results in our purely dephasing case in a cycle of length 
4At, which is identical in structure to a CP cycle (see Sec. 
IIV CI . and whose periodic repetition we term PCDD2. For 
a single qubit undergoing arbitrary decoherence, the corre- 
sponding PCDD2 protocol (constructed from a 16-pulse base 
cycle) has been shown to be the best performer in suppressing 
the effects of a quantum spin bat h 24,41 ' 43 . 

Figure [9] shows a comparison of PDD and PCDD^ proto- 
cols for i = 2, 3, for the shortest pulse separation compatible 
with the exciton qubit constraint, At = 0.1 ps. One can infer 
that, for the At and T storage values considered, PCDD^ per- 
forms better (that is, displays higher coherence maxima) than 



B. Uhrig decoupling 

We now assess the limitations of the optimal sequence pro- 
posed by UhrigiS when significant restrictions on At are in 
place. In UDD, consecutive pulses are spaced according to 




which implies, in particular, closely spaced pulses at the be- 
ginning and the end of the evolution period. Such a control 
sequence strongly suppresses the dephasing for a storage time 
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a considerable proportion of the total evolution time. This 
places a relatively large restriction on how many pulses may 
be used within a given storage period, and eventually results 
in large amounts of dephasing during the long time delays in 
which no pulses occur. 



0123456789 10 



t (ps) 



FIG. 10: | exp(— F(f))| 2 for the exciton qubit in the presence of 
UDD with a) 100 pulses, b) 50 pulses, c) 40 pulses, d) 25 pulses, 
e) 14 pulses corresponding to the best allowed sequence for the case 
of the exciton qubit. For comparison, f) shows the free evolution. 



of the order of 
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two ~ (n + 1) 



2ir : 



(38) 



where r c denotes, as before (Sec. [TTTb- the relevant bath cor- 
relation time. As mentioned, with huj c w 2 meV, this corre- 
sponds to t c « 2.06 ps. Beyond £tjddi the efficiency of UDD 
falls rapidly. From Eq. d38l ). we find that for UDD to effi- 
ciently protect the exciton qubit over T storage « 10 ps, n must 
be on the order of 100. Fig. [K)]shows the resulting UDD per- 
formance as n is decreased. It can be seen that as n < 100, 
the advantage of UDD is rapidly lost. 

For our QD system, however, the main physical limita- 
tion is on the time delay between pulses. The shortest in- 
terval between control pulses in UDD, AiJ^P , is before the 
first pulse, and after the last pulse. From Eq. fl37l i we see 
that such a sequence with n = 100 pulses over a period of 
^storage = 10 ps corresponds to Ai™ D = 2.4xl0~ 3 ps, which 
is roughly two orders of magnitude less than that allowed by 
the physical constraints for the exciton qubit in question. Even 
for a sequence consisting of n = 40 pulses only (for which 
the efficiency is already poor as shown in Fig. [10] curve c)), 
At^P D = 1.5xl0~ 2 ps, which is still an order of magnitude 
shorter than allowed. 

To respect the physical constraints, one may estimaate that 
allowed UDD sequences should have a number of pulses 
?i < 14 within the intended T^age = 10 ps. Such a se- 
quence corresponds to curve e) in Fig. [TO] It is then clear 
that any UDD sequence compatible with our physical con- 
straints is outperformed by the best allowed PDD sequence 
which would preserve a coherence close to 1 for the same time 
window (see Fig.|2]i. Fig.[10]also shows that any constrained 
UDD sequence performing like curve d) or worse would in- 
crease the dephasing compared with the free evolution, that is, 
would result in decoherence acceleration. The reason for the 
shortfalls of UDD in our setting stems from the large spread of 
the control intervals (ti — ti-i). If we impose a lower bound 
on the minimum time interval, other intervals must take up 



C. Carr-Purcell decoupling 



We now focus on analyzing more closely CPDD, which re- 
sults from the periodic repetition of a CP cycle of the formi 



At cp X2A7; cp XA7; cp . 



(39) 



This also corresponds, as noted, to PCDD2 with At cp = At 
(cf. Table IJi. Specifically, we are interested in comparing a 
PDD sequence with a CPDD having the same cycle time, T c — 
2Ai, thus At cp = At/ 2: though the corresponding pulse time 
interval may not be allowed by the physical constraints we are 
considering, this study will pave the way to be the analysis to 
be developed in the next section. 

Basically, CPDD may be viewed as a PDD protocol where 
pulses are uniformly spaced by 2Ai cp , except that the se- 
quence is displaced forward by ti = At/2, the time at which 
the first pulse is applied. As a consequence of the symmetry 
of the control propagator in Eq. d39l with respect to the cycle 
mid-point, it is well known 3 - that CPDD is a second-order pro- 
tocol as compared to standard (asymmetric) PDD, with lead- 
ing corrections of order T 3 . Using the exact representation 
established in Eq. (fT2l , we will now assess the extent to which 
CPDD improves over PDD for a purely dephasing system, and 
gain insight into asymptotic properties. 

We begin by determining the dephasing half-way be- 
tween consecutive control pulses for the case of PDD. Using 
Eq. (fT2l . we find 



T PDD \t= ln+^)At 



1 



2 J2(-l) m+1 r (mAt) 



rn — 1 



+4^^r (( m -j)A()(-ir^ 

m—2 j<m 



+2^(-l) m+ 'T 

m— 1 

+(-l)T [(n + i) 



n + — — m 1 At 



At 



(40) 
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which we may rewrite as 

pPDD [ 
n [ 



= 2^(-l) m+1 r (mAt) 



m=l 



+ 4 E E r o(( m -j) A< )(- 1 )"" 1+J 

m— 2 j<m 



k--)At 



+2]T(-i) fc+1 r c 
/c=i 

+(-l)»r [(n + i) 



I At 



(41) 



where fc = n — m + 1. Similarly, by using Eq. (fT2l . we may 
also determine the dephasing for CPDD with At cv = At/2 
and ti = (i - 1/2) At, that is, 



r PDD [(n + 1/2) At] = rf DD [(n + l/2)Ai]. (42) 

This exacf result is illustrated in Fig. QT| where we plot the 
dephasing behaviour under PDD and CPDD for the exciton 
qubit with At = 0.1 ps. As predicted by Eq. d42l , the co- 
herence in the presence of each sequence is equal at times 
t = (n+l/2)Ai. Interestingly, for t > t sat , r PDD [(n+l/2)A] 
are local maxima of coherence whereas r^ PDD [(n + l/2)At] 
are local minima, proving CPDD to be much more efficient 
than PDD provided that the time of the first pulse is allowed 
to be t\ — 0.05 ps. 




FIG. 11: Comparison of CPDD (dotted line) with At CP = 0.05 ps 
and PDD (solid line) with At = 0.1 ps for the exciton system. 



V. TOWARDS OPTIMIZED SEQUENCES IN THE 
PRESENCE OF PULSE TIMING CONSTRAINTS 

Building on the understanding gained from the comparison 
between different protocols in Sec. IIV1 we now specifically 
aim to optimize DD performance for a bosonic dephasing en- 
vironment when pulses are subject to a minimum pulse-delay 
constraint. 

The basic observation is to note that if after an initial arbi- 
trary pulse sequence, PDD is turned on at a time £pdd, then 
for t > t PDD + t sat , we have T n+1 (t + At) - T n (t) = AT^ 
(recall Eq. d25ll). This naturally suggests an interpolated DD 
approach, where an initial sequence is chosen to minimize 



r(t™ x ), whilst transforming the oscillations of coherence into 
phase with a PDD sequence to be turned on immediately af- 
terwards. Interestingly, a similar philosophy has been invoked 
to optimally merge deterministic and randomized DD meth- 
ods to enhance performance over the entire storage tima^. In 
our case, CPDD is indeed the simplest example of this inter- 
polation: as already noted, CPDD can be thought of as a PDD 
sequence applied at ipoD = At/2 + At, following a prepara- 
tory sequence consisting of a single pulse at t = At/2. 

Unfortunately, standard CPDD is not allowed in our sysyem 
due to the physical constraint: the time interval between 
pulses in the initial sequence is smaller than the minimum al- 
lowed At which characterizes the subsequent PDD sequence. 
Simply using a CPDD sequence which does not break the time 
constraint is clearly not optimal. If the smallest allowed pulse 
interval is Ai m j n , then the best CPDD sequence consists of 
periodic repetitions of a CPDD cycle with Ai cp = Ai m ; n , 
and the most efficient allowed PDD sequence is repetitions of 
X At m [ n X Ai m j n . Since CPDD cancels the terms in the Mag- 
nus expansion up to to the second order, over the first few 
repetitions it performs much better than PDD, which only can- 
cels them up to the first order. However, for longer times the 
effects due to the higher-order Magnus corrections accumu- 
late, and they turn out to do so more favorably for PDD. This 
manifests itself in a smaller Ar^ for PDD than for the best 
allowed CPDD protocol. As shown by Eqs. dB3b and (IB6I) . 
the coherence oscillations are independent of the timing of 
any pulses applied before t — t sat . Therefore, CPDD can be 
treated as a PDD sequence with At = 2At cp for t > i sat . 
This justifies defining a Ar^ for a CPDD sequence. 

Physically, what is needed is a different initial sequence that 
efficiently 'engineers' the transition of the coherence oscilla- 
tions - from the natural response frequency determined by the 
first bit-flip to the frequency of the following PPD sequence. 
To accomplish this, we propose to use CP cycles with vary- 
ing At CP . That is, we define such an interpolated sequence 
by letting the zth cycle to be characterized by a pulse de- 
lay Atf p , and begin immediately after the previous cycle at 
ti = ti-i + 4At^ p 1 . The analysis of the resulting averag- 
ing properties may be carried out by adapting the derivation 
of Ref. ,5j to the pure dephasing bosonic setting of Eq. (Q]). 
While the detail of the calculations are included in Appendix 
C, the result is that, similar to standard CP, the proposed DD 
sequence still cancels the terms in the Magnus expansion up to 
the second order. Therefore, the interpolated scheme does not 
only perform well for small t, but also quickly results in pulses 
uniformly separated by At m i n - resulting in a small Ar^, 
hence high performance for long storage times. 

The simplest way to generate a good interpolated DD se- 
quence is to apply a CP cycle with At C¥ — Ai m ; n , followed 
by periodic repetitions of one with At CF = At^/2. The 
sequence is then given by 



= At 



min s 
3At m j n 

3 Ai m i n 



-A< n 



i > 3. 



(43) 



12 



We compare this sequence with standard PDD with At = 
At m ; n in Fig. [12j upper panel. One clearly sees that the se- 
quence in Eq. (03]) is more efficient. 

By construction, the first two CP cycles in the above se- 
quence play the role of modifying the frequency of the de- 
phasing oscillations in such a way that they are brought in 
phase with the following repeated cycles. We can perform 
this process more smoothly by gradually reducing A< cp from 
At m ; n to At mm /2 over more than a control cycle. Though for 
very small Ai m j n the two cases would be equivalent, for sys- 
tems such as the exciton qubit where the time restrictions are 
relatively severe, the smoother transition sequence may de- 
couple the qubit more efficiently. Such a modified sequence 
may be implemented by applying CP cycles with decreasing 
At CP , that is, 



At' 



CP 



At m i n 
At mm 

At 



i = 1 

(i - 1)A 2 , l< i < i PDDl (44) 

/2, I > ipDD, 



where zpdd = At mm /(2A 2 ) and A 2 is an arbitrary time de- 
fined such that ipoD is an integer. The greater ipdd, the longer 
the time over which the the decreasing length cycles are ap- 
plied. Alternatively, we may describe the above sequence in 
terms of the pulse times: 



U = Ai min + (i - l)2At, 
(»-!)(» -2) 



U-i + At n 



A 2 , i < 

. ^ Atmin 
t > — : 



At„ 



A 2 
- 2. 




(45) 



t (ps) 

FIG. 12: (Color online) Comparison of PDD (dark line) with inter- 
polated DD sequences of varying CP cycles (light line), for At m i n = 
0.1 ps. Top: Sequence given by Eq. d43b - Bottom: Sequence given 
by Eq. i44\ , with A2 = 0.01 ps. Inset: Zoom over the initial part of 
the time window with timings of the pulse sequences explicitly indi- 
cated (diamonds for PDD and squares for the modified sequences). 



The above modified sequence is compared to PDD with At = 
Atmin in Fig. Q~2] lower panel. As before, we see that DD 
with varying CP cycles outperform PDD. Furthermore, it can 
be seen that there is an improvement over the more abrupt 
sequence described by Eq. d43l . 



In Sec. IIV A 21 we showed that a constrained CPDD se- 
quence could outperform PDD over time scales of the order 
of 10 ps (see top panel of Fig. [9] for PCDD 2 ) even though, 
for longer times, the smaller Ar^ for PDD would eventually 
make it more efficient than CPDD. In Fig. [13] we further com- 
pare CPDD with the interpolated sequence given in Eq. d44b . 
The latter is found to be slightly more efficient than the best 
allowed CPDD sequence over short time scales. Furthermore, 
because of the smaller Ar^, as time progresses it will also 
outperform CPDD asymptotically. A main advantage of the 
sequence given in Eq. (|44] |. however, is that it not only leads to 
higher maxima than CPDD, but also, after the first few pulses, 
to a much smaller coherence oscillation amplitude. This re- 
flects the fact that the oscillation period has been tuned to the 
minimum allowed time interval At m ; n . In this respect, the per- 
formance of the sequence in Eq. ( PPfl i is more robust against 
the precise readout times, or, equivalently, a readout offset er- 
ror relative to the coherence maxima would not significantly 
affect the coherence recovered using this sequence. 



VI. CONCLUSION AND OUTLOOK 

We have investigated the ability of DD to inhibit decoher- 
ence of a single qubit coupled to a purely dephasing bosonic 
environment, by comparing the performance of low-level pe- 
riodic DD schemes based on uniform pulse separations to 
higher-level non-uniform DD schemes. 

For arbitrary spectral density functions characterizing dif- 
ferent dephasing environments, we have derived an exact rep- 
resentation of the controlled dynamics available for instanta- 
neous pulses, and exact relationships for the asymptotic dy- 
namics in the case of PDD. Building on these results, we have 
shown that a main weakness of PDD is due to the oscilla- 
tion of coherence following the first bit-flip being out of phase 
with the rest of the sequence. This has naturally suggested the 
application of a suitably engineered preparatory sequence as 
a strategy to enhance DD efficiency, by bringing the coher- 
ence oscillations into phase with a subsequent PDD sequence. 
The resulting 'interpolated' DD protocols are found to be es- 
pecially efficient for physical systems where the mimimum 
time interval between control pulses is strongly constrained. 
For such systems, DD protocols like concatenated or Uhrig 
DD, which are designed to achieve peak performance when 
the asymptotic regime of arbitrarily small pulse separations is 
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FIG. 13: (Color online) Comparison of the best allowed CPDD se- 
quence (light line) with a sequence of CP cycles of decreasing length 
(dark line) as described in Eq. d44t with A2 = 0.01 ps and for 
At m in = 0.1 ps. Inset: Zoom over the initial part of the time window 
with timings of the pulse sequences explicitly indicated (squares for 
the modified sequence and crosses for CPDD, respectively). 



fully accessible, tend to largely lose their advantages. 

For the excitonic dephasing environment of interest, in par- 
ticular, we have shown how a sequence of Carr-Purcell cycles 
with suitably chosen (analytically generated) time delays pro- 
vides a very efficient DD protocol for realistic QD parame- 
ters and qubit storage times. Our process of constructing a 
DD sequence under which the coherence oscillates asymp- 
totically with the minimum period allowed by the physical 
contraints offers, as by-product, the advantage of a signif- 
icantly smaller coherence oscillation amplitude, relative to 
constrained PCDD or UDD sequences. This makes the pro- 
posed interpolated sequence more robust against readout. 

While our analytically-designed interpolated DD proto- 
col might be compelling in its simplicity, identifying DD 
schemes that are guaranteed to yield optimal performance 
subject to non-trivial timing constraints appears as an interest- 
ing control-theoretic problem for further investigation. Revis- 
iting the local numerical optimization approach recently pro- 
posed in Ref. [If] in a constrained minimization perspective 
might offer a concrete starting point in this respect. Like- 
wise, the investigation of dynamical error-control schemes 
based on bounded-strength 'Eulerian' DD^S, along with the 
recently proposed extension to decoherence-protected quan- 
tum gate s 47 ' 48 , might prove especially fruitful for exciton 
qubits, in view of the reduced control overheads associated 
with purely dephasing environments. Lastly, an interesting 
general question is to what extent exact representations of the 
controlled coherence dynamics in terms of the uncontrolled 
one may exist for an arbitrary purely dephasing error model, 
allowing, for instance, exact insight into the long-time con- 
trolled dynamics to be gained. 
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APPENDIX A: DERIVATION OF Ar o 



For PDD, t n = nAt, using this expression and recasting 
the sums in Eq. ( l22l in terms of k = n — j, we can write 



Ar 



PDD 



(-l) n r ((n + l)Ai) - 3r (nAi)(-l) n 



4£r (*At)(-i) 



(Al) 



k=l 



Using Eq. (0 for n = and extending the sum to include 
k = 0, Eq. ( IA1 b becomes 

/>oo />oo 

Ar PDD = _ 4 / v ^duj + (-l) n 6rj{(j)cos(wnAt)duj 
Jo Jo 

poo 

-(-1)" / 2ij(oj)cos(uj(n + l)At)dw 
Jo 

poo n—l 

+8 / ij{uj)^2cos(ujkAt){-l) k duj. 



(A2) 



By using the relationship 



A, lk , ( 1 (-l)»cos(^±Ix) 
Y.i- 1 ) c os(kx) = --+ hr—t, (A3) 



fc=i 



2 cos 



we can rewrite the above equation as 



Ar 



PDD 



T}(w) {(-l) n+1 cos[(n + l)ujnAt] 



+ (-l) n+i cos(nujnAt) 
+ 2 (-D"cos(yA tM )1 



X cos 



(2n+ 1 



AtuA 



f Atuj\ 



\ 2 
(_1)« cos(^tlAtw) 




(A4) 
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This finally can be rearranged as 



Ar 



rj(uj) sin 2 



ujAt 



(-l) n cos [u (n + i) At] 



du. (A5) 



We now take the limit of n — > oo, and note that 

(-l) n cos(w(n+±) At) 1 



lim 

n — >oo 



COS 



27T 



— ; (5 (uj - (Wres + 2Zcj res )) 



At 



(A6) 



If we assume that r/(aj res ) S> f7(w res + 2/w r es) for Z > 0, which 
is true for sufficiently small At (that is, w res 3> w c )> we can 
neglect contributions from Z > and define 



AIV = lim Ar 



PDD 



(A7) 



8dujS(u! — uj Tes )r](Lij)L0 ms sin 2 



uAt 



8r](u> res )u) ri 



APPENDIX B: LONG-TIME LIMIT OF AT™ D (t) 

By using Eq. ( fT~2T > and straightforward manipulations, we 
can separate Ar™ D (t) from Eq. ( f24b into two parts, 



Ar^© = Ar^ + Ar-(t), 



with 



Ar™(t) = (-1)" [r (*„ + t) - r (t n+1 + *) 

and the second term 



(Bl) 



(B2) 



ATI 1 = 2(-l)' 



r (i„+i) + 2^(-i) J r (i„+i - 1,) 



independent of i. By using 

r (t) = 2 / - cos(w<)]dw, 

we can rewrite Ar™(£) as 

/•OO 

Ar™(t) = -4(-l)" / , ? Hsm 
Jo 

{[(' 



(B3) 



(B4) 



sm • 



+ -)Ab + t 



u^djj. (B5) 



The last term in the integrand above is fast oscillating for large 
n, so we will have that for n > n sat 



where e can be made arbitrarily small. 

Let us now consider AT™. By using Eq. (1B41 i. the relation 
Eq. ( IA3b and some tedious but straightforward manipulations, 
we can rewrite Eq. dB3b as 



Ar 



TI 



-(-1)"4 / n(aj) cos[w(n + l)Ai]du 
Jo 

4 / 77 lj - — - Vat-^ '-du. (B7) 

7o 2 cos A*- 



Again, the integrand in the first term of the above equation 
is fast oscillating for large n, while the second term tends to 
AEoo for n — ► 00 (see Eq. dA61 l. We can then write that for 

n > n sat , 



Ar! 1 - Ar^i <e. 



(B8) 



By combining Eqs. dB6l ) and dB8l ), we finally obtain that for 

any t and n > n sat , 



Ar™ D (t) « Ar^. 



(B9) 



We note that in the case of supraohmic environment, by 
using that r (oo) = lim,,-,^ r (t„+i) = 2 f™ri(u))du is 
finite, and using Eqs. (IB2l i and ( IB7l i. we can recast the condi- 
tions Eqs. (IB6b and dB8l ) as 



(BIO) 



|r (t>* Mt )-r (oo)| <e, 



where t sat = n s . dt At. This emphasizes that condition 
applies for times at which the natural evolution saturates to 
its long-term behavior. 



APPENDIX C: AVERAGING PROPERTIES OF 
INTERPOLATED DD SCHEME 



We begin by casting the QD Hamiltonian Eq. ([]]) (with 
a = 1/2) in the following form: 



Hi = a z ® B z + cr ® B , 



(CI) 



where B z and Bq are operators acting on the phonon bath, and 
do and a z denote the identity and the Pauli matrix acting on 
the exciton qubit, respectively. This allows us to express the 
evolution in the presence of the i-th CP cycle by the propaga- 
tor 



,CP\ 



U^{AAtY) = U f {AtY)XU f {2AtY)XU f {At 

where Uf(t) = exp(— ti?i) represents free evolution for a 
time t. If we define 



AT^(t) \<e, 



(B6) 



Ho 



-a z ® B z + cr <E> B = XHiX, 



(C2) 



15 



we can write the entire sequence propagator as a Magnus se- 



ries expansion^, 



(C3) 



where H (t) = U^(t)HU ctI i(t) is the time-dependent (piece- 
wise constant, for instantaneous pulses) effective Hamilto- 
nian that describes the evolution under the control propagator 
U ca \(t) resulting from the applied pulses^^- 



for which, in the limit of sufficiently fast control, we can only 
consider the first two lowest-order terms in A*£. Specifically 
(in units where h = 1): 



Ax 



t 

i I dtiH(ti), 
ft r ti 



(C4) 



A 2 = -- / dh dtiMhlHfa) 



For the sequence of different CP cycles described in Sec. 
IVl A\ is proportional to the identity operator, and hence does 
not contribute to dephasing. This is a simple consequence 
of the qubit spending equal amounts of time in each of the 
(C5) computational basis states. More interestingly, we find 



An 



t«+4Ai cp 



r tf+4At cp 



f 1 _ _ c 

dti dt 2 [H(t 1 ),H(t 2 )]=y2 dh I dt 2 [H(ti),H(jk2)\ 

Jo i=1 Ju Jo 

V / dti dt 2 + dh dt 2 + dti dt2)[H(t 1 ),H(t 2 )}\ 

LV -'* -'t i +At cp JO Ju+3At ce JO ' J 



E 



ti+At cp i-1 

diiJ^Aff [H U H 2 ] 



3 = 1 



JO Jti+3At c 
t l +3At cp / i-1 

dt!^2Atf + At?)[H 2 ,H 1 ] 



ti+At cp 



ti+4At c 



i,+3A( cp 



dti^2Atf [Hi, Ha]} 
j'=i 

£ { AtP ^ 2Atf [Hi, ff a ] + 2AtP C |] 2A«f + Atf p ) [ff 2 , ^] + Aif ^ 2A<f [Hi, F 2 ]} 



= o, 



i=i 



again up to irrelevant pure-bath terms. This confirms the second-order cancellation claimed in the main text. 
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